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We present a search for the decay B — > r v T using 467.8 x 10 6 BB pairs collected at the T(4S) 
resonance with the BABAR detector at the SLAC PEP-II B-Factory. We select a sample of events 
with one completely reconstructed B~ in an hadronic decay mode (£>~ — > Z)'*' X~ and B~ — ► 
J/ipX^). We examine the rest of the event to search for a B + — > t + v decay. We identify the 
t + lepton in the following modes: r + — > e + ^ e 77 T , r + — > ^v^Vt, t + — s> tv + V t and r + — > pv T . 
We find an excess of events with respect to expected background, which excludes the null signal 
hypothesis at the level of 3.3 a and can be converted to a branching fraction central value of 
B{B + -> t+u t ) = (1.80tojI(stat.) ± 0.26(syst.)) x 10~ 4 . 
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INTRODUCTION 



The study of the purely leptonic decay is of particular interest as a test of the Standard Model (SM) and a search 
for physics beyond the SM. It is sensitive to the product of the B meson decay constant fs, and the absolute value 
of the Cabibbo-Kobayashi-Maskawa matrix element \V u b\ [l|. In the SM the branching fraction is given by: 



B(B+ 



f B \Vu b \ gn 



1 



TB+, 



(1) 



where Gf is the Fermi constant, t b + is the B + lifetime, and m B and m T are the B + meson and r lepton masses. 

The process is sensitive to possible extensions of the SM. For instance, in two-Higgs doublet models and in 
minimal supersymmetric extensions of the SM it can be mediated by a charged Higgs boson. A branching fraction 
measurement can therefore also be used to constrain the parameter space of extensions of the SM. In a previously 
published analysis, based on a tagging technique using hadronic B decays that is similar to that used in this paper and 
a smaller data set, the BABAR collaboration measured B(B + t + v t ) = (1.8+°;g(stat.)±0.4±0.2(syst.)) x 10" 4 Q, and 
using tagging based on reconstruction of semileptonic B decays B(B + — > t + v t ) = (1.7± 0.8(stat.) ± 0.2(syst.)) x 10~ 4 
The Belle collaboration measured, with a similar tagging technique used in this analysis, the branching fraction to 



beB{B+ -> t+v t ) = (1.79t^(stat.)+°;5i(syst.))xl0 
of semileptonic B decays B(B + 



[H , and using a tagging algorithm based on the reconstruction 



v T ) = (1.54±g;f 7 (stat.)+°f (syst.)) x 



10" 4 @. 



THE BABAR DETECTOR AND DATASET 



The data used in this analysis were collected with the BABAR detector at the PEP-II storage ring. The sample 
corresponds to an integrated luminosity of 426 fb -1 at the T(4S) resonance (on-resonance) and 44.51b" 1 taken at 
40MeV below the BB production threshold (off- resonance) , which is used to study background from e + e~ — > f f 
(/ = u, d, s, c, r) continuum events. The on-resonance sample contains (467.8 ± 5.1) x 10 6 BB decays. The detector 
is described in detail elsewhere 0- Charged particle trajectories are measured in the tracking system composed of a 
five- layer silicon vertex detector and a 40-layer drift chamber (DCH), operating in a 1.5 T solenoidal magnetic field. 
A Cherenkov detector is used for charged ir-K discrimination, a Csl calorimeter (EMC) for photon and electron 
identification, and the flux return of the solenoid, which consists of layers of iron interspersed with resistive plate 
chambers or limited streamer tubes, for muon and neutral hadron identification. 

In order to estimate signal selection efficiencies and to study physics backgrounds, we use a Monte Carlo (MC) 
simulation based on GEANT4 8]. In MC simulated signal events one B + meson decays as B + — » t + v t and the other 
decays in any final state. The BB and continuum MC samples are equivalent to approximatively three times and 1.5 
times, respectively, the accumulated data sample. Beam-related background and detector noise are taken from data 
and overlaid on the simulated events. 



SIGNAL SELECTION 



We reconstruct an exclusive decay of one of the B mesons in the event (which we refer to as the tag B) and 
examine the rest of the event for the experimental signature of B + — > t + v t (charged-conjugate modes are implied 
throughout the paper). We consider the most abundant r decay modes r + — > e + vv 1 t + — > \v~vv, t + — > 7r + z^, 
t + — > p + v, totaling approximatively 70% of all r decays. The signal region in data is kept blind until the end of the 
analysis chain when we extract the signal yield. We reconstruct the tag B candidate in the set of hadronic decays 
B~ — > M°X~, where M° denotes a Z?M° or a J/ip, and X~ denotes a system of hadrons with total charge — 1 
composed of nxir , n2K ± 7 n37r°, 714K® where n\ + ri2 < 5, < 2, 713 and m < 2. We reconstruct the D° as 
D" K-'K + ,K-T: + 'K a ,K-T: + n-7r + ,K^T: ,K°T: + TT- 1 K°TT + T:~T: ,K + K-,T: + 'K'. We reconstruct the D*° meson as 
jj*a _^ D°ir° , D°j, and the J/ip meson as J/ip — > e + e~,fi + /i~. The kinematic consistency of the tag B candidates 
is checked with the beam energy-substituted mass toes = ^/s/4 — p' B and the energy difference AE = Eb — \fs/2, 
where yfs is the total energy in the T(4S) center of mass system and ps and Eb denote respectively the momentum 
and the energy of the tag B candidate in the center of mass frame. The resolution on AE is measured to be 
a AE = 10— 35MeV, depending on the decay mode; we require \AE\ < 3a/±fi. Events with a candidate tag B arise 
from two possible classes with different toes distributions. Signal events with a correctly reconstructed tag B and 
the other B decaying as B rv , and background events from T(AS) —> B + B~ with a correctly reconstructed tag 
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B are characterized by an toes distribution peaked at the B mass. The other class of events consists of continuum 
background, e + e~ — s- qq (q = u, , s, } and e+e~ — > t + t~, and combinatorial background, T(4S') — > B°B° or B + B~ in 
which the tag £? is misreconstructed; this class of events has a broad toes distribution that can be modeled by means 
of a phenomenological threshold function (ARGUS function) Q . 

If multiple tag B candidates are reconstructed we select that with the lowest value of \AE\. The purity V of each 
reconstructed B decay mode is estimated as the ratio of the number of peaking events with tties > 5.27 GeV to the 
total number of events in the same range. We consider only events with the tag B reconstructed in decay modes with 
V > 0.1. The yield in data is determined by means of an extended unbinned maximum likelihood fit to the toes 
distribution, as shown in figure[T] We use as probability density function (PDF) for the combinatorial and continuum 
background an ARGUS function, while for the correctly reconstructed tag B component we use as PDF a Gaussian 
function with an exponential tail (Crystal Ball function) 10]. Combinatorial and continuum backgrounds in any 
discriminating variable are estimated from a sideband in toes (5.209 GeV < toes < 5.260 GeV) and extrapolated into 
the signal region (toes > 5.270 GeV) using the results of a fit to an ARGUS function. The peaking B + B~ background 
is determined from B + B~ MC, after subtraction of the combinatorial component to avoid double counting by means 
of a similar fit. 
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FIG. 1: Fit to the ?71es distribution in data. Dots are data, the blue curve represents the fitted combinatorial and continuum 
background. 



After the reconstruction of the tag B, we apply a set of selection criteria on the rest of the event. We require 
the presence of only one well reconstructed charged track (signal track), with charge opposite to that of the tag B. 
The r lepton is reconstructed in one of four decay modes: r+ — > e + vD, r+ — > p + v&, r+ — > ir + v, r+ — > p + v. We 
separate the event sample in four categories using particle identification criteria applied to the signal track. The 
t + — > p + v sample is obtained by associating the signal track with a ir° reconstructed from a pair of neutral clusters 
with invariant mass between 115 MeV/c 2 and 155 MeV/c 2 . In order to remove the e + e~ — > t + t~ background we 
impose r mode dependent requirements, preserving 90% of the B + — > t + v t signal, on the ratio between the 2 nd and 
the th Fox- Wolfram moments (R2) [ll[ calculated using all the charged tracks and neutral clusters of the event. 

In order to reject the continuum and combinatorial background we use discriminating variables constructed 
from the kinematics of the tag B candidate. The first variable is the momentum in the CM frame (p* M ) °f 
£)(*)o or jjr^ candidate reconstructed from the decay products of the tag B. The second variable is the abso- 
lute value of the thrust 12| (|Tb|) of the tag B. The third variable is the cosine of the angle between the thrust 
of the tag B and the thrust of the rest of the event (cos^ts). We combine p* M , \Tg\ and cos 6tb in a likeli- 
hood ratio L c = Ls(p* M , \Tb\, cos 9 T b)/(Ls(p* m , \T b \, cos6 T b) + L B {p* M , \T B \, cosOtb)), where the signal (S) and 
background (B) likelihoods are obtained from the products of the PDFs of the three discriminating variables: 
Ls(p* M ,\f B \,cos6 TB ) = Ps(p* m )Ps(\Tb\)Ps{cos9tb) and L B (p* M , \f B \, cos9 TB ) = P B (p* M )P B (\f B \)P B (cos 8 TB ). 
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The PDFs for the signal modes are obtained from the signal MC, whereas the PDFs for backgrounds are obtained 
from the toes sideband in data. 

In order to further reject the background from correctly reconstructed tag B events, we impose a requirement on 
center of mass momentum of the signal track for the r+ — > e + vv , r+ — > p + vv and t + — > tt + p modes. For the 
t + — > p + v mode we combine in a likelihood ratio (Lp) the following variables: the invariant mass of the signal track 
and the 7r°, the total momentum in the CM frame of the pair \p* p \, the momentum in the CM frame of the tt°, and 
the missing mass of the event. The PDFs used in the likelihood ratio for the signal and background are determined 
from signal and B + B~ MC, respectively. 

The most discriminating variable is -E ex tra, defined as the sum of the energies of the neutral clusters not associated 
with the tag B or with the signal tt° from the t + — > p + v mode, and passing a minimum energy requirement (60 MeV). 
Signal events tend to peak at low -E C xtra, background events, which contain additional sources of neutral clusters, tend 
to be distributed at higher values. 

We optimize the selection requirements, including those on the purity V of the tag B and the minimum energy 
of the neutral clusters, aiming at the lowest expected uncertainty in the branching fraction fit. In order to estimate 
the uncertainty, which includes the statistical and the largest systematic uncertainties, we run 1000 toy experiments 
extracted from the background and signal expected shapes for a set of possible selection requirements, assuming a 
branching fraction of 1.4 x 10~ 4 jl3||. 

The signal selection requirements are summarized in table |TJ The -E ox tra distribution with all the selection require- 
ments applied is shown in figure [2j 



Variable r + — > e + vv r + — > p + vv t + — > -k + v t + — i p + v 

purity > 10% 

cluster energy ( MeV) 60 

R2 < 0.57 < 0.56 < 0.56 < 0.51 

L c > 0.2 > > 0.3 > 0.45 

p* trk (GeV/c) <2.1 <2 >1.4 

L P >0.8 



TABLE I: Optimized signal selection criteria for each r mode. 



BRANCHING FRACTION MEASUREMENT PROCEDURE AND RESULTS 

We use an extended unbinned maximum likelihood fit to extract the B + — > t + v branching fraction. The likelihood 
function for the N k candidates reconstructed in one of the four r decay modes k is 

C k = e -(«-.»+»».*) JJ J ns.kVm.k) + n b , k V b k (E itk ) \ (2) 
i=i ^ ' 

where n s .k is the signal yield, n^^ is the background yield, 2?^ is the E extlll value of the i th event, Vf. is the probability 
density function of signal events, and V\ is the probability density function of background events. The background 
yields in each decay mode are permitted to float independently of each other in the fit, while the signal yields are 
constrained to a single branching ratio via the relation: 

n s , k = N B g x e k x B (3) 

where N B g = (4.678 ± 0.051) x 10 8 is the number of BB pairs in the data sample, e k is the r decay mode dependent 
reconstruction efficiency, and B is the B + — > t + v branching fraction. The parameters N B -g and e k are fixed in the fit 
while B is left floating. The reconstruction efficiencies e kl which include the r branching fractions, are obtained from 
MC simulation of the signal. Since the tag B reconstruction efficiency is included in e k and is estimated from the 
signal MC, we apply a correction factor i?data/MC = 0.926 ± 0.010 to take into account data/MC differences, taking 
the ratio of the peaking component of the toes distribution of the hadronic tag B in data and in MC simulation 
events. 

We use histograms with a bin width of 60 MeV to represent the PDFs 7-^ and V\ for signal and background, 
respectively. The signal PDF is obtained from a high statistics signal MC simulation sample, corrected for data/MC 
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FIG. 2: -Bextra distribution in data (dots with error bars) with all selection requirements applied and fit results overlaid. 
The hatched histogram is the background, the red dashed component is the signal. Plot (a) shows all r decay modes fitted 
simultaneously. Lower plots show the projection of the simultaneous fit result on the four analyzed r decay modes: (b) 
t + — > e + vu , (c) t + — > n+vv , (d) t + — > Tv + f , (e) r + — ► p + v . 



disagreement. Since a data sample of suitable statistics with exactly the same final states as our signal channels is not 
available, we use a sample of fully reconstructed events where in addition to the reconstructed tag B, a second B is 
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Decay Mode 


ex 10 


Branching Fraction (xlO J Sif 


mificance a 


— 7 1 C^~/V/V 
T + — S> [l^VV 


2.73 
2.92 


1 oq+0.89 
1 - zo -0.80 


0.5 
1.6 


T + — ¥ 7T + f 


1.55 


4.oll:i 


3.3 


T + — ¥ p + U 


0.85 




2.6 


combined 


8.05 


1.80^ 


3.6 



TABLE II: Reconstruction efficiency e, measured branching fractions and statistical significance obtained from the fit with all 
the modes separately and constrained to the same branching fraction. The r decay mode branching fractions are included in 
the efficiencies. 



reconstructed in an hadronic or a semileptonic decay mode, using charged tracks and neutral clusters not assigned 
to the tag B. In order to estimate the correction to the signal PDF, we compare the distribution of -E e xtra in this 
double tags sample from experimental data and MC simulation. The distributions are normalized to the same area 
and the comparison is shown in figure [3] We extract the correction function taking the ratio of the two distributions 
and fitting it with a second order polynomial. 

We take the PDF of the combinatorial background from the Toes sideband. The contribution of this component 
in the signal region is obtained by fitting the toes distribution after the selection has been applied. The shape of 
the peaking background is taken from B + B~ MC in the signal region, after the intrinsic combinatoric background 
has been subtracted by a fit to toes, to avoid double counting. The two background components are added together 
in a single background PDF. We finally apply a smoothing procedure on the total background shape, excluding 
the first bins. We estimate the branching fraction minimizing — 21n£, where C — n^. =1 jC,fe , and C% is defined in 
equation OThe projections of the fit results are shown in figure [2j We observe a significant excess of events with 
respect to the expected backgrounds and measure a branching fraction B(B + -> T +v T ) = (1.80i°;^) x 1CT 4 , where 
the uncertainty is statistical. We evaluate the significance of the observed signal, including only statistical uncertainty, 
as S — v/21n(jC s +(, / '£&) , where £ s +{, and C s denotes the obtained maximum likelihood value and the likelihood value 
assuming background only. We find S = 3.6a. Table ITT1 summarizes the results from the fit. 

SYSTEMATICS 

The dominant source of systematic uncertainty is the background PDF, due the finite statistics of the B + B~ MC 
simulated sample, used to estimate the B + B~ background PDF, and of the toes data sideband, used to estimate the 
combinatorial and continuum backgrounds. In order to estimate the systematic uncertainty we repeat the fit of the 
branching fraction with 1000 variations of the background PDF, varying each bin within the statistical error, and 
assign 12% as systematic uncertainty. 

The systematic uncertainty due to the signal E'extra distribution correction function obtained from data/MC com- 
parisons using control samples is obtained by varying the parameters of the second order polynomial within their 
uncertainty and repeating the fit to the B + —> t + v t branching fraction. We observe a 1.7% variation that we take as 
the systematic uncertainty on the signal shape. 

Uncertainty in the differences between data and MC in the tracking and neutral reconstruction efficiencies reflects 
in the uncertainty in the central value of the branching fraction. The difference of the tracking efficiency is estimated 
with a control sample of high momentum tracks from e + e~ — > t + t~ events to be 0.5% per track. Since there is only 
one signal track candidate in all four r decay modes in the B + — > t + v t signal, we use this value as the uncertainty 
due to tracking efficiency. We accept events with one extra low px charged track. Comparing the multiplicity of low 
Pt charged tracks from the double tags sample in data and in MC, we estimate the systematic uncertainty to be 1.3%. 
Adding in quadrature the two uncertainties we estimate the the systematic error to be 1.4%. 

Other systematic uncertainties on the efficiency stem from the finite signal MC statistics (0.8%), the uncertainty 
in the tag B efficiency correction (5.0%), the electron identification (2.6%), and muon identification (4.7%). The 
systematic uncertainties are summarized in Table IIIII The total systematic uncertainty is obtained by combining all 
sources in quadrature. 
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Hadronic DT, E [GeV] Hybrid DT, E [GeV] 



FIG. 3: -Eextra distribution for double tags. The second B is reconstructed in hadronic decays (left plot) or semileptonic decays 
(right plot). Points are data, histograms are MC simulation. 



Source of systematics BF uncertainty (%) 



B counting 


0.5 


Tag B efficiency 


5.0 


Background PDF 


12 


Signal PDF 


1.7 


MC statistics 


0.8 


Electron identification 


2.6 


Muon identification 


4.7 


Kaon identification 


0.4 


Tracking 


1.4 


Total 


14 



TABLE III: Contributions to systematic uncertainty on the branching fraction. 



CONCLUSIONS 

In summary, we have measured the branching fraction of the decay B + — » t + v t using a tagging algorithm based on 
the reconstruction of hadronic B decays using a data sample containg 467.8 x 10 6 BB pairs collected with the BABAR 
detector at the PEP-II £?-Factory. We measure the branching fraction to be B(B+ ->■ t+v t ) = (1.80t^(stat.) ± 
0.26(syst.)) x 10~ 4 , excluding the null hypothesis at the level of 3.6 standard deviations using statistical uncertainties 
only, and at the level of 3.3 standard deviations including the systematic uncertainties. This result supersedes 
our previous result using the same technique Combiningthis result with the other BABAR measurement of 
B(B + — s- t + v t ) derived from a statistically independent sample [4|, we obtain a single BABAR result B(B + — > t + i> t ) = 
(1.76 ± 0.49) x 10~ 4 , where the uncertainty includes both statistical and systematic uncertainties. 
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